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ABSTRACT 

Wc observed the accreting millisecond pulsar SAX J1808. 4-3658 with Gemini-South 
in g' and i' bands, nearly simultaneous with XMM-Ncwton observations. A clear 
periodic flux modulation on the system's orbital period is present, consistent with 
the varying aspect of the donor star's heated face. We model the contributions of 
a disk and donor star to these optical bands. To produce the observed modulation 
amplitudes, we conclude that the donor must be irradiated by an external flux 2 orders 
of magnitude greater than provided by the measured X-ray luminosity. A possible 
explanation for this irradiation is that the radio pulsar mechanism becomes active 
during the quiescent state as suggested by Burderi et al., with relativistic particles 
heating the donor's day-side face. Our modelling constrains the binary inclination to 
be 36-67°. We obtain estimates for the pulsar mass of > 2.2M Q (although this limit 
is sensitive to the source's distance), consistent with the accelerated NS cooling in 
this system indicated by X-ray observations. We also estimate the donor mass to be 
O.O7-O.11M0, providing further indications that the system underwent non-standard 
binary evolution to reach its current state. 

Key words: binaries : stars: individual (SAX J1808. 4-3658) - stars: pulsars - stars: 
neutron - X-rays: binaries 



1 INTRODUCTION 

The low-mass X-ray binary (LMXB) SAX J1808.4- 
3658 (hereafter J1808) is the first detected accretion- 
powered millisecond pulsar dChakrabartv fc Morgan! 1 19981 : 
IWiinands fc van der Klisl Il998h . Such objects are thought 
to be evolutionary intermediaries between LMXBs and ra- 
dio millisecond pulsars (MSPs): they are still in an LMXB 
phase with ongoing mass transfer but have had time for 
mass transfer to spin up the accreting neutron star (NS) 
to millisecond spin periods. If and when accretion halts 
in these s ystems, they are expected to turn-on as ra- 
dio M SPs (|Alpar et al.l 1 19821 : iRadhakrishnan fc Srinivasanl 
Il982t ). Since J1808's discovery, seven more accretion- 
powered millisecond pulsars have been discovered, but J1808 
stands out as the best observed member of this class. 

J1808 has offer ed us more than several surp rises. Its 
NS's mass function (|Chakrabartv fc Morgan! 1 19981 ) implies 
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that the donor co mpanion (which is almost cert ainly a low- 
mass brown dwarf lBildsten fc Chakrabartvl200lj ) has a mass 
> O.O4M0 (assuming a fiducial NS mass of 1.4 Mq), too 
large to be explained b y standard bi nary evolution scenar- 
ios (see, e.g., figure 5 of[Dcloyc 2008|). The system's orbital 
period, Port., appears to be evolving an order of magni- 
tude faster than bina ry theory predicts (|di Salvo et al.ll2008l ; 
lHartman et al. I l2008l . although, see their discussion of pos- 
sible explanations for this). X-ray observations of the sys- 
tem during its quiescent phases are only able to place upper 
limits on the NS's thermal component, indicating the NS 
cools extremely rap idly as compared to most other NSs (see 
I Heinke et al.l 120071 and references therein) . Perhaps most 
surprisingly, observations indirectly suggest J 1808 transi- 
tions between the L MXB and radio MSP state during X- 
ray quiescent phases ([Homer et ai1l200ll ; iBurderi et al 1l2003l : 
ICampana et al.ll2004h . 

The disk in J 1808 is thermally unstable and undergoes 
outbursts lasting approximately a month roughly every 2 
years. During the inter-outburst quiescent phase (when the 
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disk is cooler and dim and may not extend inward close 
to the NS), the system exhibits roughly sinusoidal variabil- 
ity in the optical with the same 2 hr Pp r b de tected in X- 
rays l|Homer et al.ll200ll ; ICampana et al] 12004 ) . The phas- 
ing of optical maxima corresponds to whe n the donor is di- 
rectly behind the NS (|Homer et al.|[200ll ). indicating that 
the donor's face is heated by flux originating near the NS. 
Typically this flux would be due to the X-ray radiation pro- 
duced by the accretion flow through the inner accretion disk 
and onto the NS surface. However, the amplitude of the opti- 
cal modulation requires significantly greater X-ray flux than 
is observed during quiescence (albeit the X-ray observations 
are non-simultaneous with the optical). This led to spec- 
ulation that th e NS commences radio pulsar activity dur - 
ing quiescence (|Burderi et alj 120031 ; ICampana et al.l l2004h , 
providing the necessary amount of flux to heat the donor's 
day-side face. 

To investigate J1808 in both optical and X-ray wave- 
lengths nearly simultaneously, we have obtained imaging op- 
tical data (g' and i') from the Gemini Observatory separated 
by 6.5 hours from an XMM observation of SAX J1808 in qui- 
escence in 2007. Our science goals for these observations in- 
cluded further constraints on the thermal component of the 
X-ray emission and its X-ray variability in quiescence, si- 
multaneous determination of the optical modulation to con- 
strain definitively its origin, and determining if optical data 
allows constraints on either the quiescent disk's structure or 
on the binary's parameters. The X-ray observations are de - 
scribed in detail in a companion paper l|Heinke et al.ll2008T ). 
Here we discuss the Gemini observations. In 32 we describe 
the observations and data reduction. In Sj3] we describe the 
light curve analysis within the framework of a theoretical 
model. We then discuss the constraints on this model's pa- 
rameters obtained from our data in i|3.4l In iJH and O we 
discuss our results and summarize. The Appendix includes 
details on how our data and further system details (in par- 
ticular J1808's distance) constrain our model's parameters. 



2 OBSERVATIONS 

We observed 1808 on March 8 and March 10, 2007 with 
the Gemini Mu lti-Object Spectrograph-South (GMOS-S; 
iHook et al.ll2004l ) in imaging mode. The March 10 observa- 
tion (program GS-2007A-Q-8, UT 06:47 to 09:55) was nearly 
simultaneous with the XMM-Newton observation (ObsID 
0400230501), which started at 16:24 UT, and continued for 
16 hours without evidence for strong variability. We took a 
time series of 44 g' exposures of 185 s each, with two 230 
s i! exposures at either end to obtain colour information. 
Atmospheric seeing was generally good, ranging from 0.65" 
to 0.98" FWHM. Landolt standard star fields PG1047+003 
and SA104 were also observed during the same night. The 
March 8 observation, using an identical program, suffered 
from poorer seeing (FWHM ranging from 1.0 to 1.4"), but 
still provided data adequate for our purposes, and thus we 
analyse it as well. 




Figure 1. Finding charts for J1808, average of the seven g' 
GMOS-S frames with the best (0.65") seeing. A nearby (0.5") 
star can be barely distinguished to the SE. 



2.1 Optical photometry 

We used the calibrated Gemini data products, which are 
processed with Gemini-specific IRAF □ tasks to remove the 
bias, and flat-field and mosaic the images. We performed 
hotometry on g and i' reference frames using DAOPHOT 
Stetson 1987) point-spread-function (PSF) fitting photom- 
etry, calibrated using the standard star fields. 

Compari s on of our images with the finding chart of 
iHomer et alj (|200l[ ) reveals two stars, separated by 0.5", 
that are consistent with J1808's indicated location (see Fig- 
ure [I]). The brighter of these two is variable on a 2-hour 
period, and we confidently identify it with J1808. By asso- 
ciating ten uncrowded unsatur ated nearby stars w ith stars 
in the USNO B1.0 catalogue jMonet et alj|2003l ). we find 
a position of a=18:08:27.63, <5=-36:58:43.37 (J2000), with 
uncertainties of 0.2" in each coordinate (accounting for 
the uncertainty in the transformation to the USNO B1.0 
frame). This i s con sistent with the VLA-derived position of 
iRupen et all d2002l). and with th e optical position recently 
derived by Hartman et alj (|2008r i. The optical position de- 
rived by iGiles et alj ( 19991 ) in outburst, 1.7" away, is also 
marginally consistent when the absolute errors in the GSC 
(up to 1.6" in the southern h emisphere. iTaff et al.|[l990l ). 
used as the reference frame bv lGiles et al] ~( 19991 ) . are con- 
sidered. 

The fainter star (g'=22.35±0.04, i'=20.68±0.02) is lo- 
cated only 0.5" from J1808 to the SE, at a=18:08:27.67, 
#=-36:58:43.66. This star has been identified in previous un- 
published images of the J1808 field (D. Chakrabarty 2007, 
priv. comm.), and may affect psf-fitting or aperture photom- 
etry of J 1808 in quiescence taken in poor seeing. 

We then used the differential photometry program ISIS 



1 IRAF (Image Reduction and Analysis Facility) is distributed 
by the National Optical Astronomy Observatories, which are op- 
erated by Association of Universities for Research in Astron- 
omy, Inc., under cooperative agreement with the National Science 
Foundation. 
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Figure 2. Calibrated g' light curves for J1808, from March 10 
(crosses) and March 8 (triangles). The fit shown i s to a sinusoid 
with orbital period fixed to the value of 7249.1569 s dPapitto et al.l 
l2005ft . 



(|Alard fc LuptorJll99St I Alardl l2000h for e ach fi lter, follow- 
ing the basic method of iMocheiska et ail (|2000l ) as follows. 
We select the best-seeing reference frame for each filter, and 
transform the remaining frames onto the reference coordi- 
nate systems, aligning the frames with bright stars. We pro- 
duced a reference frame from the seven best-seeing g' frames. 
This frame is convolved with a spatially varying (3 degrees 
of freedom) kernel to match each frame's PSF, and sub- 
tracted from the remaining frames. Only variable stars (and 
stars which are saturated) should remain on the subtracted 
images, as stars of constant brightness cancel out. Finally, 
profile photometry is extracted from the subtracted images, 
barycenter corrected to TDB using the Perl extension As- 
tro::Time::HJD to correct the observation times to helio- 
centric Julian dates. We performed aperture photometry for 
the standard stars, and used this photometry to convert our 
PSF-fitting photometry, and our differential ISIS photome- 
try, into calibrated g' and i' magnitudes. 



3 ANALYSIS OF THE OPTICAL LIGHT 
CURVES 

3.1 Data Processing 

The observed optical light curves show significant variability 
on time scales of order the orbital period as well as flickering 
on time scales shorter than P or b. Having only obtained four 
i' values during each night's observations precludes averag- 
ing the i' data between nights; the intrinsic variability of the 
source further precludes a simple combination of i' data into 
a single light curve. We proceed by analysing each night's 
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Figure 3. The final form of the March 08 and March 10 optical 
data used in our model light curve fits. Two complete orbits are 
shown for clarity. 



light curves separately. Later we'll consider what constraints 
on J1808's system parameters can be derived by considering 
the two nights' data in tandem. 

In order to compare against our model light curves, 
we sorted o ur data by orb i tal p hase, 4> using the orbital 
ephemeris of iPapitto et all (|2005l Fl The March 10 g data 
exhibits a clear, systematic brightening beyond the second 
observed minimum (i.e., past ~ 9.1 hrs in figure [2]|, as well 
as large amplitude flickering. This aperiodic variability is 
likely due to stochastic processes in the disk (e.g, hot-spot 
variability) that are not captured by our models described 
below. We do not have a sufficient number of observations 
to average over the stochastic variability. So, in order to 
compare our models with a observational time-frame repre- 
senting as much as possible almost constant conditions in 
the system, we excluded from our analysis the March 10 
(/-data beyond the second flux minima. This keeps the g' 
points obtained during the same orbit as the March 10 i' 
data. The remaining g data we phase sorted, binned (using 
a bin width A(f> = 0.05), and then averaged by bin to re- 
duce the flickering's impact on the fits. The data resulting 
from all these procedures, and used in our model light curve 
fitting below, are shown in Figure [3] 



3.2 Theoretical Light Curve Modelling 

We utilize the program ELC (|Orosz fc Hauschildt|[2OO0l ) to 
calculate theoretical model light curves. To model the opti- 
cal data, we consider contributions from both an accretion 
disk and a donor that is irradiated by an energy flux originat- 
ing at the location of the accreting NS. The ELC program 
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includes the effects of limb darkening (implicitly via use of 
PHOEN1X0 model atmospheres) and disk/donor occupa- 
tions. Our model has seven independent input parameters. 

The first two are the masses of the accreting NS, Mi, 
and the donor, Mi. Once Mi, M2 are specified, the binary's 
inclination relative to our l ine of sight, ibin is determine d 
via the NS's mass function (|Chakrabartv fc Morgan! 1 19981 ). 
We determine the binary's orbital separation from Mi, Ms, 
and Porb- From P or b and the requirement that the donor fills 
its Roche lobe, the donor's radius, R2, is determined once 
M2 is specified. This also fixes the donor's unirradiated ef- 
fective temperature, T e s, which we obtain from unpublished 
low-mass brown dwarf models of Deloye & Taam (the T e g 
obtained from thes e models are very similar to those of, e.g., 
iBaraffe et al.lll998h . The resulting T e s, which are a function 
of M 2 alone, range between fa 2000-3000 K. 

To account for the periodic variability, the donor's op- 
tical contribution must be phase dependent. The ELC pro- 
gram allows modelling this source of variability by irradiat- 
ing the donor with a point source of flux at the accretor's 
location. We consider the impact of such a point source with 
total luminosity, Li rr , between 10 33 and 10 34 ' 6 erg s _1 . We 
fix the donor's albedo in ELC to the value minimizing the re- 
quired Pirr for a given light curve amplitude 0. If the donor's 
albedo differ from this, a commensurate increase in Li rr will 
be required to produce the same results. 

To model the disk in J1808, we use four parameters: the 
disk's inner and outer radii, ri n and r ou t, the temperature of 
the disk at n B , Tin, and an exponent, £, that characterizes 
the disk's radial temperature profile, T(r) tx Ti n (r/ri n )^ . 
The structure of a quiescent accretion disk, in particular 
is not expected to be that of the standard steady-state 
disk model (e.g.. iDiibus et al]|200ll ). This expectation has 
been verified by observations of other binaries in quiescence, 
which find, i n particular, £- values ranging from £ < — 1 to 
£ « (e.g.. ISkidmore et al.1 12000|; IVrielmann et all [20021 ; 



IShahbaz et all 120031 ; iBaptista fc Bortolettol |2004 steady 
state disks are expected to have —0.75 < £ < —0.5). Given 
this expectation, we treat each of these four parameters as 
independent quantities. 

The values of Li rr used here are substantially greater 
than t he Lx ~ 10 32 ergs s -1 observed from 1808 in quies- 
cence |Heinke et al.l [2008h . As seen below, such high val- 
ues of Li rr are required to fit our data, making it clear 
that the X-ray emission in the system insufficient to drive 
the obs erved optical variab i lity. As has already b een sug- 
gested (|Burderi et al.1 120031 ; ICampana et alj|20t)4 >. the NS 
itself could provide the required irradiating flux in the form 
of a pulsar-wind that turns-on during the quiescent phase 
of the disk outburst cycle. However, we will refrain from 
detailed attempts to interpret the source or quality of the 
required irradiating energy. Instead, we view our results as 
a guide to the amount of flux required to impinge on the 
donor in order to produce the observed amplitude of vari- 
ability. We will make a few further comments on this point 
in the discussion. 

We also attempted to fit the optical data using an 
irradiated-donor only model. For this case, instead of fix- 
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ing the donor's unirradiated T e g at the values predicted by 
brown dwarf evolution models, we allowed T e ff to be a free 
parameter ranging between 4500 and 7500 K (chosen to pro- 
vide sufficient g' flux at light curve minimum given J1808's 
distance as discussed below). This provided an input pa- 
rameter set of Mi, M2, T e ff, and Li rr . This model was moti- 
vated by the possibility that horizontal fluid currents in the 
donor's atmosphere may efficiently advect the irradiating en- 
ergy to the donor's night side, elevating the surface tempera- 
ture t here (see, e.g. jBurkert et aLll2005 ; lDobbs-Dixon fc Linl 
l200&t ). However, this model can not simultaneously repro- 
duce the light curves' morphology, amplitude, and i! — g' 
colours, and did not produce acceptable fits. 

We used ELC to calculate model light curves at discrete 
values of six of our input parameters: Mi = 1.4, 1.6, 1.8, 
3.0 M ; M 2 = 0.05, 0.06, 0.07, 0.13 M Q ; r in = 0.001, 
0.01, 0.03, 0.06, 0.1, and 0.2 P L ,i (where R L>1 is the Roche 
lobe radius of the accretor and is the internal unit used by 
ELC ); r out = 0.25, 0.4, 0.5, 0.6, 0.7, 0.8, and 0.9 R L ,i; Tin = 
5000, 5500, 6000, 10000 K; and £ = 0.0, -0.1, -0.2, -0.3, 
-0.5, and -0.75. 

For each set of input parameters, V — (Mi, M2, rin, 
r ut, Ti n , £), and each night's data, we first minimized the 
X 2 contribution of the g' data, Xg', against Li IT . We did 
not calculate the i' data's contribution to \ 2 , X%, simul- 
taneously with the x 2 ,' -determination because ELC ignores 
colour information during simultaneous fits. Thus, once we 
determined the L 1TT that minimized x%', we calculated the 
corresponding i! light curve and calculated x\> retaining the 
i' — g information. We finally calculated x 2 = Xg' + Xi'- 

We found that to obtain "acceptable" fits (that is, a 
reduced x 2 ~ 1), we needed to systematically increase the 
data's error bars by 0.0164 and 0.0129 mags for the March 
08 and March 10 data, respectively. Both of these values are 
typically less than the la errors shown in Figures [5] The 
error bars shown in Figure [3] reflect this increase. 

To validate that our spacing of the input parameters 
is fine enough to effectively probe the x 2 hyper-surface, we 
performed additional fits using several test grids of input 
parameters with spacings at least twice as fine as above. The 
best fits obtained in these test calculations were no better 
than on our standard grid. 



3.3 Possible Influence of the Accretion Stream 
Hot Spot 

The ELC code does not model contributions from the hot 
spot where the accretion stream impacts the disk. We argue 
here that the hot spot's potential contribution can not ex- 
plain the observed light curves. However, it may contribute 
to the disk's overall flux at an important level; if this were 
the case, then the disk parameters we derive from our fits 
would be suspect. Given that we are unable to meaningfully 
constrain the parameters of our disk model, as discussed in 
H3.4I this shortcoming does not impact our central results. 

The main argument against a modulation of the hot 
spot's contribution explaining the observed light curve is 
the orbital phase of flux maximum. We can compare the 
time of the ascending node phase, T aBC , from RXTE timing 
with our Gemini optical data. If the sinusoidal modulation 
represents the heated companion-with light maximum cor- 
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responding to donor's superior conjunction-then the T a3C 
from the Gemini data is (TDB) MJD=54169.2784(3). We 
can compare this with the pr edicted T asc advance d from the 
combined orbital solution of IPapitto et al.l (|2005T l. which is 
(TDB) MJD=54169.280(6), wit h which it is consisten t. Us- 
ing the more recent ephemeris of lHartman et all (|2008h gives 
a predicted T asc of (TDB) MJD=54169.28000(3), differing 
from our T asc by 138±47 seconds, or 1.9 ± 0.6% of the or- 
bital period. Either ephemeris confirms that light maximum 
occ urs when the do nor i s behi nd the NS. 

lHartman et~afl (|2008l ). which is (TDB) 
MJD=54169. 28000(3). These times differ by 138 ± 47 
seconds, or 2 ± 1% of the orbital period, confirming light 
maximum occurs when the donor is behind the NS. 

On the other hand, the hot spot maximum should lead 
the donor's inferior conj unction by some de gree. Extrapo- 
lating the calculations of iFlanneni (|l975T l or lLubow fc Shul 
|l975l ) to the expected mass ratio range of J1808 predicts 
that, for a disk extending only out to the stream's circu- 
larization radius, the hot spot maximum should lead the 
inferior conjunction by a phase of 0.1-0.13. I.e., at a phase 
0.37-0.40 later than the observed maximum. The circular- 
ization radius represents the minimum distance the hot spot 
can be located from the NS. The disk likely extends beyond 
this radius, leading to a hot spot location that leads the 
donor's inferior conjunction by an even smaller margin. This 
strongly supports the origin of the modulation as due to the 
heated secondary, rather than the hot Spot. 



3.4 Parameter Constraints from Light Curve Fits 

3.4-1 Model Constraints from, Source Distance Estimates 

Before estimating input parameter values from the above 
X 2 calculations, we utilized J1808's distance estimates to 
only consider theoretical m odels whose (/-band flu x es fal l 
within appropriate ranges. iGallowav fc Cummind (120061 ) 
constrained J1808's distance, d, using several different meth- 
ods. They modelled the type I X-ray bursts that occurred 
during the 2002 October outbursts, comparing X-ray burst 
recurrence times, bolometric fluence, and the ratios of in- 
tegrated persistent flux to burst fluence. From their best 
fitting burst models, they estimate that d = 3.1-3.8 kpc. 
The observed X-ray bursts all showed evidence for radius 
expansion; for a pure helium NS atmosphere (i.e., con- 
sistent with the X-ray bur sts being He triggered events, 
IGallowav fc Cumminej I2006TI . this constrains d = 3.6 kpc 
(this estimate decreases if there is any H present in the atmo- 
sphere) . Finally, based on the fluenc es and recurrence times 
of the system's long-term outbursts, IGallowav fc Cummind 
(|2006t ) estimate a lower limit of d > 3.4 kpc by assuming 
the minimum time-averaged mass transfer rate onto the NS 
is set by the binary's gravity-wave driven M in the conser- 
vative limit. 

Based on these estimates, we take as a fiducial con- 
straint d^_3 ; 5Jcp_c_^_a_3%jm£ertamty (i.e., the d = 3.4-3.6 
kpc of lGallowav fc Cumm ing 2006). For comparison, we also 
carried out our analysis assuming an uncertainty of ±10% 
on d. We then used these d estimates to convert the observed 
gr'-values at light curve minima into a (/-band luminosity at 
4> = 0.0, L g i^ =0 . For t his we took A a i — 0.85 based on 
Nh = 1-29 x 10 cm - 2 (jDickev fc Lockman|[l990l '); the cor- 



responding Aii = 0.47. The resulting L g i ranges for the 
3% (10%) uncertainties on d are 2.13-2.42(1.82-2.74) x 10 30 
and 2.47 - 2.80(2.11 - 3.16) x 10 30 erg s" 1 for the March 08 
and March 10 data, respectively. 



3.4-2 Data and L g i ^^Q-Limit Influences on Input 
Parameter Constraints: General Trends 

There are several general trends in how the input parame- 
ters influence the fits. First, the disk parameters (n n , r out , 
Ti n , and £) are all completely unconstrained by the obser- 
vations. Since each of these parameters can influence over- 
all disk brightness and colour, the fact they are not indi- 
vidually constrained is not too surprising. Second, x 2 > is 
mainly a function of i D m. The p'-data on its own (see dis- 
cussion below) rules out ibin > 70° at the 3-<r level (i.e., 
A\ 2 — \ 2 > — Xmin > 9-0 f° r one f fee parameter). Third, 
the best-fitting L 1IT increases with L g i and, in general, 
Z/irr > 10 34 ergs s _1 is required to produce acceptable fits. 

The L g i it j, =0 limits indirectly constrain both ibin and 
Mi. The relative contribution of the disk vs. the donor- 
which impacts x 2 primarily through a model's i'-band 
amplitude-impacts both L 9 / ^ =0 and Mi. Larger L 9 / ^ =0 - 
values have larger (phase independent) disk contributions, 
resulting in a lower i'-band amplitude. Larger Mi-values 
have (for a given ibin) larger and hotter donors, which tend 
to increase i'-band amplitudes. These two effects in combi- 
nation lead to minimum x 2 " va lues at increasing L 9 / ^ =0 as 
the Mi-value is increased. Within the L g / =o -range derived 
from the distance constraints, the minimum x 2 - vames gen- 
erally decrease with increasing Mi while the minima-x 2 at 
fixed Mi are increasing functions of L g /^ = o (see upper panel 
of Fig. IA4|) . Thus the lower L o /0 =o -limit indirectly favours 
more massive NSs. On the other hand, the high-L g / ^=o limit 
indirectly disfavours lower values of ibin since models with 
brighter L g i^ =0 in general require larger M2. 



3-4-3 Input Parameter Constraints 

Taking our fiducial d constraint, we are able to produce ac- 
ceptable fits to each night's data: the best-fit to the March 
08 data had a x 2 = 38.70 for 37 degrees of freedom (d.o.f), 
while the best-fit to the March 10 data had a x 2 = 23.94 
(for 17 d.o.f.). Taken individually, each night's observations 
place constraints on ibin and Li rr at the 2a level. Both 
nights require ibin > 30.0 at a 2a confidence. For March 08, 
log(L irr ) = 34.14-34.25 erg s" 1 , while the March 10 limits 
log(L irr ) = 34.06-34.20 erg s -1 . The March 08 i' data, since 
it includes points near both light curve maximum and min- 
imum, constrains our models more strongly then the March 
10 data. Because of this, the March 08 observation, on its 
own, provides limits on J1808's component masses, requiring 
Mi > 1.8M and M 2 > O.O7M (where these limits are on 
the two parameters individually). The preference for larger 
Mi results from the need to achieve the observed i'-band 
amplitude given the lower limit on L g /^ =0 required by the 
d-estimates (see discussion in i]3.4.2|l . 

While the March 10 data set does not by itself provide 
strong constraints on J1808's binary parameters, fits to it do 
still favour larger Mi values. Thus we also considered both 
nights observations in tandem in order to determine if more 
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Figure 4. The model light curves that, in combination, provide 
the best fit to both night's data at fixed Mi, M2 within our fidu- 
cial Lgi ^—Q limits. The solid lines show the theoretical models. 
The symbols show the data (same as in Fig. [3j used to calculate 

x 2 . 



stringent constraints on Mi and M2 could be developed. To 
do so, for each pair of (MijMO-vahies, we considered all 
input-parameter sets Vk (producing a fit with Xk against 
March 08 data) and Ve (producing a fit with \t against 
the March 10 data) sharing the specified (Mi,M2)-values. 
As the source distance is treated as a free parameter by 
ELC when determining x 2 , we needed to be careful to only 
consider (Vk, T^-pairs whose resulting best-fit distance es- 
timates, f^ELC , were mutually consistent between the two 
nights. We considered a (Vk, Ve) to have mutually consis- 
tent d E LC if \d,ELC ,i - rfELC ,j| < Ad 3(T , where Ad 3(T is the 
maximum change in model distance producing a A\ 2 < 9-0 
in a fit's quality. We determined Ads a by averaging over the 
required change in cJelc for « 60 of our best-fitting models 
between both nights. Finally, for each (Vk, Ve) with mutu- 
ally consistent distances, we calculated x\,e = Xk + Xe- 

For our fiducial L g i ^ =0 limits, the best fitting summed 
Xl,e = 63.21 (for 57 d'.o.f.) at Mi = 3.0M Q and M 2 = 
O.IM©. The light curves of the two models combining to pro- 
duce this best fit are shown in Figure [4] In Figure [5] we plot 
the difference between the minimum x\ e at each (Mi,M2) 
and this best-fit x\ ^-value as grey-scale values keyed to con- 
fidence limits. One can see from Fig [5] that our results con- 
strain Mi > 2.2M and 0.07 < M 2 < 0.11M Q at the 2a 
level; at the same confidence level, 36 < ibm < 67°. Thus, 
our optical data provides evidence for J 1808 harbouring a 
massive NS, consistent with expectations from the apparent 
rapid cooling the NS in this system undergoes. 

The derived lower limit on Mi does depend on the lower 
Lg'.tp^o limit, which depends on d. Thus, to explore how 
sensitive our parameter estimates are to the lower L g i 
limit, we also determined limits on Mi and M 2 assuming a 
10% and 20% uncertainty on d. The results for the case of 



A X 2 




M 2 (M ) 

Figure 5. Parameter constraints derived by combining fits to 
both the March 08 and 10 data. The grey squares show the dif- 
ference, Ax| 1, between the minimum value of x\ e a * eacn (Mb 
M 2 ) and the global x 2 ^-minimum. The greyscale values indi- 
cate the uncertainty levels associated with this Ax| t for a sin- 
gle independent parameter (i.e., one and two a corresponds to a 
Ax| i = 1.0 and 4.0, etc.); thus various regions provide limits on 
Mi and M2 taken independently. The solid white and black lines 
provide approximate contours of the 2<x and 3<x levels to guide 
the eye. For this plot, the uncertainty on J1808's distance, d, was 
our fiducial 3%. 

a 10% uncertainty are shown in Figure [6] Under these as- 
sumptions, we constrain Mi > 1.8Mq, M2 > 0.06Mq, and 
32 < ibin < 74° . For the case of a 20% d-uncertainty, we can 
not place any limits on Mi or M2. Upcoming Gemini time 
series observations of J 1808 will provide better colour infor- 
mation, and may provide stronger evidence for a massive 
NS. 



4 DISCUSSION 

4.1 Comparison to Prior Work 

There have been several prior observations o f SAX J1808 in 
the optical during quiescence. iHomer et alj (|200ll ) reported 
on observations in both white light and BVR filters. Each 
band's light curve was modulated on the system's P or b and 
they estimated the variabi lity in the l^-band to h ave a semi- 
amplitude of ~ 0.08 mag. ICampana et al. I (|2004l) report on 
/-band photometry of J1808, again observing light curve 
modulation on the s ystem's P or b. They also re-evaluated the 
IHomer et al.l l|200ll ) data, reporting a U-band modulation 
semi-amplitude of 0.13 ± 0.06 mag. 

For our best-fitting models, the (/-band modulation 
semi-amplitude for either night is ~ 0.3 mag. Since these 
prior observations did not resolve the fainter star near 
J1808's position, the fact that we find a larger modulation 
amplitude is not surprising. Quantitatively, this faint star 
produces « 66 — 80% the p'-band flux of J1808 at light curve 
minimum. If our observations had not resolved this star, we 
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Ax 2 




M 2 (M ) 

Figure 6. Same as Fig. [5] but for the case of a 10% uncertainty 
on d. 



I. To proceed along these lines, we determined lower limits 
on Li rr as a function of M\. These limits are derived from 
the March 08 data due to its larger (/-amplitude by taking 
the minimum March 08 Li rl within the set of all (Vk, Vi) 
with specified M\ and a xi,e within 3-ct of the minimum. 

For our fiducial 3% d-uncertainties, this minimum 
Li rr ~ 1.5 x 10 34 erg s _1 , almost independent of Mi. With 
v = 400.975 Hz and v = 5.6 ± 2 x 1(T 16 Hz s _1 , this leads 
to I > 1.7 ± 0.6 x 10 45 g cm 2 (with the quoted error on / 
only including the v uncertainty contribution). The central 
value of this / limit is certainly consistent with our finding 
that J1808 harbours a more massive NS and, for a range of 
nuclear eq uations of state, wo uld require M\> 1.4Mq (see 
Figure 4 of lWorlev et al.|[2008h . However, the large v uncer- 
tainty prevents our Li rr limits from excluding much of the 
7- Mi parameter space and even Mi « 1.2M© are consistent 
with the lower end of the allowed /-range. We note in pass- 
ing that if the efficiency at which the pulsar's particle wind 
energy is thermalized in the donor's atmosphere were 50%, 
then the lower limit on I, including the v uncertainty, would 
rule out Mi < 1.8 Mq at least fo r the r ange of equations of 
state considered in IWorlev et al.l |200ct ). 



would have found (/ -band semi-ampl i tudes of « 0.15 — 0.19, 
consistent with the Campan a et all (|2004T ) analysis of the 
V-b and data oflHomer et all (|200ll ). 

IWang et al. (|200ll ) used optical observations taken dur- 
ing J1808's 1998 outburst to constrain the system's ibi n and 
Ay. For an assumed distance of 3 kpc, their 90% confidence 
interval on i is 20-60° , with a best-fitting value of 37° . Our 
2a ibin limits (36-67°) are consistent with these results. 



4.2 Estimates of the Neutron Star's Moment of 
Inertia 

Since our observations were able to resolve the faint inter- 
loping star nearby J1808, it is not surprising that we need 
a larger L[ lr to explain the larger optical modulation am- 
plitudes we find. Quantitatively, in our good fitting models 
Lj r r is greater by a factor of ~ 2.5 — 3 relative to that found 
bv lCampana et all |2004 ) . As with these prior observations, 
our required Li lr is 1-2 orders of magnitude greater than 
the observed quiesc ent X-ray luminosity ~ 8 x 10 31 erg s _1 
|Heinke et al.ll2008l ). We also checked whether the amount 
of flux required to irradiate the donor could be provided by 
the remnant disk in the system and found that the donor is 
only able to reprocess < 1% of the disk's luminosity. Thus, 
this scenario falls 1-2 orders of magnitude short of being able 
to explain the modulation amplitudes. We conclude that the 
only source that plausibly can power the optical modulation 
is the spin-down ene rgy of the central NS. As suggested by 
iBurderi et all (120031 ). this energy is likely tapped through 
a pulsar-wind that turns-on during quiescence (when the 
mass-transfer rate onto the NS is reduced enough to allow 
the accretion disk to be truncated outside of the NS's light- 
cylinder) . 

The luminosity provided by the NS's spin-down is E = 
An 2 Iuu, where / is the NS's moment of inertia and v is 
its spin- frequency. iHartman et al.l (|2008h have recently mea- 
sured J1808's long-term v = -5.6 ±2 x 10" 16 Hz s" 1 . Thus, 
the requirement L m < E allows lower limits to be placed on 



4.3 Neutron Star Cooling Rate and Mass 
Estimates 

Our Mi estimates can also be co nnected with the NS coo l- 
ing rate observed in J1808 (e.g.. iHeinke et al.1 120071 , 120081 ). 
Observations indicate that NSs exhibit a range of cooling 
rates across both the isolated and accreting sub-populations 
that are set by the specific neut rino emission processes oc - 
curring in the NS core (see, e.g.. lYakovlev fc Pethickll2004l ). 
Generally speaking, the most efficient neutrino production 
process available to lower mass NSs is the modified Urea 
process, which result s in a "standard" cooling rate (see, e.g., 
lYakovlev et~aT1l200ll ). At the higher core densities found in 
higher mass NSs, several direct Urea processes-depending 
on the constituent nature of matter at high density (i.e., 
nucleons, hyperons, or "free" quarks)- are allowed that pro- 
vide orders of magnitude faster cooling (l Yakovlev fc Pethickl 
|2004| ) than the standard rate. 

The NS mass at which the transition from standard to 
enhanced cooling occurs and how broad the range of masses 
over which it occurs depends on the equation of state of 
matter at supra-nuclear densities and other physics (such as 
the occurrence of transitions to various superfluid states and 
in-medium modifications to the neutrino emiss ion processes 
lYakovlev fc Pethickll2004l ; lBlaschke et al.ll2004h . Most young 
isolated NSs appear to cool at rates close to the standard 
one. On the assumption that most of the NSs in the observed 
cooling sample have masses consistent with other NS popu- 
lations, we can infer that those N Ss cooling near the stan- 
dard rate have masses ~ 1.35M0 |Thorsett fc Chakrabartvl 
119991 ). Within the isolated NS population there are sev- 
eral clear cases (the Vela pulsar and PSR J0205+64) of 
objects with cooling rates intermediate between the stan- 
dard and the fully-developed enhanced rates. Effective mod- 
elling of all the isolated systems strongly indicates that 
the transition to enhanced cooling has to occur over some 
mass range and these intermediat e cooling-rate systems have 
masses in this transi tion interval (|Yakovlev fc Pethickll200l ; 
iBlaschke et al.ll2004l ). 
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Amongst the set of accreting NS systems, there is a 
similar spread in observed cooling rates. Again, most sys- 
tems cool at rates consistent with or slightly faster than the 
standard rate. However, J1808 is one of two systems host- 
ing an accreting NS whose cooling rate is consistent with 
the fully-developed enhanced rate. Our constraints limiting 
Mi > 1.8Mq are consistent with the picture of high-mass 
NSs accessing the enhanced cooling mode. We will postpone 
a detailed analysis of how J1808's NS cooling rate combined 
with Mi estimates can constrain matter's equation of state 
at supra-nuclear densities for a fut ure paper. 

We should note that recently (|Leahv et al.ll2008T ) have 
provided estimates for Mi in J1808 based on modelling the 
pulsed X-ray light curve observed during its 1998 outburst. 
Across a set of various models, they find best-fitting Mi 
values < I.IMq. Such low Mi values are inconsistent with 
J1808's enhanced cooling rate given the propensity of data 
supporting "normal" NS masses w 1.35M (see above ) . The 
3a contours quoted for some of the lLeahv et al.l (|200St ) mod- 
els do include Mi values up to ~ I.6M0, which could be 
consistent with J1808's rapid cooling. At this same confi- 
dence level, o ur fiducial Mi esti mates are barely consistent 
with those of iLeahv et~ai1 (|2008l ). 



4.4 Prior Binary Evolution and Current Donor 
State 

Give n its NS mass functi on, P OT h, and presence of hydro- 
gen jCampana et al.ll2004r ). J1808 is most likely the descen- 
dent of a NS-low mass MS binary whose evolution paral- 
leled that of the related WD accretin g cataclysmic vari- 
ables (see, e.g., iKolb fc Baraffd Il999l ). The initial evolu- 
tion of J1808's LMXB phase was towards shorter P or b as 
its then main-sequence (MS) companion was able to main- 
tain thermal equilibrium and contract under mass loss. As 
M2 decreased, its Kelvin-Helmholtz time increased, eventu- 
ally becoming longer than its mass-loss time and the donor 
began to expand under mass loss, driving the binary to 
longer P or b. Standard evolutionary models find that the re- 
sulting P or b minimum should occur near 70 minutes (while 



Kolb fc Baraffel 

(jPolitano et al.lll99Sl : 



the o bserved value is closer to 80 minutes 
19991) at which point M 2 « 0.065M o 
Kolb fc Baraffd 1 19991 ) . By the time s ystems have evolved 
back out to P orb = 2 hr, M 2 = 0.02 M m JPolitano et al]l998t 



As already noted by l|Bildsten fc Chakrabartvl 
200 ll ). the minimum M% allowed by the NS's mass function 
(O.O43M for Mi = 1.4M ) is greater than this expected 
M2-value. Our 2a limit on the donor's mass, M 2 > O.O7M0, 
provides independent evidence for M2 values significantly 
higher than standard theory predicts for J1808. 

At J1808's P orb , if the donor truly has M 2 > O.O7M , 
it would have to have a greater entropy content than the 
M 2 « 0.06 - O.O9M donors at minimum P orb a; 70 - 80 
minutes. This possible discrepancy could be explained in 
at least two ways: (a) J1808's minimum P or b was larger 
than predicted and the MS companion at this point still had 
sufficient entropy to account for the donor's current state; 
and/or (b) the donor has been heated during the system's 
post Porb-minimum evolution. The former could occur if ei- 
ther the donor's KH time was longer than expected (e.g., 
due to strong external irradiation from the accretion flow) 
or its mass-loss time was shorter than expected. This could 



be achieved via an additional angular momentum sink be- 
sides gravity-wave radiation losses or to non-conservative 
mass-transfer ef fects. Certain models for t idal heating of 
the donor (e.g.. lApplegate fc Shahaml ll994i) could provide 
a mechanism for heating the donor during the post Porb- 
minimum evolution. We plan to perform detailed modelling 
of J1808's prior binary evolution to quantify which of these 
effects can plausibly explain the system's current state. 



5 SUMMARY 

We undertook optical observations of SAX J1808 in quies- 
cence in g' and i' with Gemini South that were taken nearly 
coincident with XMM-Newton X-ray observations. The ob- 
servations were taken on two separate nights with coverage 
of slightly more than one of J1808's ~ 2 hr orbits each night. 
We detected optical modulations in both night's data on the 
system's P or b in both bands. Our observations resolved a 
nearby interloping star that was unresol ved in the observa- 
tions used in prior lig ht curve analyses (|Homer et al 1 l200ll : 
ICampana et al.ll2004 ). The optical modulations in our data 
have a larger amplitude than those measured in these prior 
efforts, and this amplitude difference is consistent with prior 
observations being contaminated by the nearby faint star. 
Our data also exhibits significant inter-orbit variability. 

We modelled the g and i' light curves obtained using 
the light curve modelling program ELC . Our model's free 
parameters were the NS and donor masses (Mi, M2), pa- 
rameters describing the disk (ri n , r out , Tin, and £), and the 
luminosity irradiating the donor, L irr , that powers the op- 
tical modulation. The binary inclination, ibm, and orbital 
separation are then determined from the measured NS mass 
function and system P or b- To fit our data, we constructed 
a grid over (Mi, M 2 , Tj n , ri n , r ou t, £) and for each vertex 
of this grid, optimized Lirr, performing the most extensive 
modelling of J1808's optical behaviour to date. 

We derived estimates for our input parameters from 
the subset of our models whose g' luminosity was consis- 
tent with the optical data given J1808's distance, d, con- 
straints. At the 2<j confidence level, we constrain ibm = 36- 
67°. We also constrain Mi and M 2 , but these limits are 
sensitive to the uncertainties in d (see the Appendix). Fo r 
our fiducial 3% d-uncertainty (G alloway fc Cummindl2006l ) . 
we constrain Mi > 2.2M and 0.07 < M 2 < O.11M . For 
a 10% d-uncertainty, these limits expand to Mi > 1.8 and 
M2 > O.O6M . We treat all the disk parameters as indepen- 
dent (since a quiescent disk is not in a steady state), and 
this freedom prevented us from deriving any constraints on 
the disk's structure. 

To power the optical modulation, L 1IT = 1.15 — 1.78 x 
10 34 erg s _1 . The almost coincident X-ray observation de- 
termined the X-ray luminosity of the source at th is epoch 
to be 7.9 ± 0.7 x 10 31 erg s" 1 (|Heinke et alj l2008h . far be- 
low that required to power the op tical variability. How ever, 
the measured NS spin-down rate (|Hartman et al.ll2008l ) can 
reason ably provide the re quired Lirr- As suggested previ- 
ously (|Burderi et al.ll2003l ). this could indicate that radio 
pulsar activity switches on in J 1808 during its quiescent pe- 
riods and the resulting pulsar particle wind provides the 
necessary flux to illuminate the donor's day-time side. If the 
donor is less than 100% efficient in converting this particle 
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flux to thermal radiation, then the required Li rr increases 
above our quoted range. An efficiency of 50% would require 
a sufficiently large NS moment of inertia that Mi < I.8M0 
would be ruled out for many reasonable nuclear equations 
of state. 

Our estimates of Mi > 1.8Mq are co nsistent with 
J180 8's very low thermal X-ray luminosity ()Heinke et al.l 
120081 ). Such massive NSs are expected to cool rapidly due 
to direct Urea neutrino emission processes that can not oc- 
cur in lower mass NSs. Our estimates that M2 > 0.07Mq 
also provide independent support (in addition to the NS's 
mass function) that the donor in this system has a signifi- 
cantly higher entropy than expected from standard binary 
evolution models. This would indicate that any of several 
processes that could increase the system's minimum P or b or 
heat the donor have been operative in the system. 
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APPENDIX A: FURTHER DETAILS ON HOW 
THE DATA AND DISTANCE CONSTRAINTS 
DETERMINE INPUT PARAMETERS 

Here we provide some of the background details concerning 
how the data and L g i >( f, =0 limits contribute to constraints 
on input parameters. Figure IA1I shows, for the March 08 
data, how the minimum of x 2 at specified Mi varies as a 
function of L 9 / ^ =0 . The solid coloured lines show the total 
X 2 , while the coloured dashed and dotted lines show the g' 
and i' data's contributions to the total x 2 (both the g' and 
total values are decreased by 10.0 to reduce the vertical scale 
of the plot). Due to the fact that we minimize x 2 ' against 
Li r r first and then determine xi afterwards, the g' data's 
amplitude is relatively well fit across all models and there is 
little variation in xt> with L g >^ = o or Mi. 

All the differences in x 2 with L g i =o therefore are pro- 



duced by the i 
different values of L 



9' ,0= 

data. For each Mi, xi 



has a minimum at 



9',0=O> 



with the Lgi^-Q at these minima 
1 , larger 



increasing with Mi . Above L g , : ^ =0 > 2 x 10 30 ergs s 
Mi values produce better fits to the i' data. The vertical 
dashed (dotted) lines show the March 08 limits on L g i ^ =0 
for our 3% (10%) uncertainties on d. We see that increasing 
the lower limit of L g i^—o removes preferentially the best fit- 
ting models at lower Mi. This is the origin of the distance 
limit producing a preference for a massive Mi. 

Before going into more detail about why larger Mi val- 
ues lead to better fits in the L g i !( j, =0 range, we'll discuss how 
the data sets the ibm constraints. In Figure IA2I we show 
the minimum x 2 versus it,i n , again at fixed values of Mi 
for the March 10 data. Colours and line styles in this plot 
have the same meaning as in Fig. IA1I and the data only in- 



clude those systems satisfying the L 



g>, 0=0 



constraints. The 



g data in this case, since we have a full light curve, on its 
own constrains the orbital inclination to a reasonable degree. 
The i' data, on the other hand for the most part provides 



a weaker trend versus 



'bij 



except at ihh 



< 30°. The in- 



creasing Xi> at the lower ibin range is actually contributed 
to by the L g i^ =Q upper- limits. Higher L g / ^ =0 values tend 
to have better fits with more massive Mi values. Removing 
the brighter systems from the sample tends to remove the 
best fitting models at lower ibin ■ Thus in combination, the 
g light curve morphology and the upper-limits on L g i^ =0 
favour ibin in the range 30 < ibin < 70°. Unlike the case for 
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Figure Al. The evolution of minimum x 2 a t fixed Mi versus 
L g i ^—q for the March 08 data. Solid lines show the overall x 2 - 
minima as a function of L g i ^— , while the dashed and dotted 
coloured lines show the g 1 and i' data's' contributions to the over- 
all x 2 - The vertical dashed and dotted lines provide the L g i l( i=o 
limits on the March 08 data assuming a 3% and 10% uncertainty 
on d. 
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Figure A2. The variation of minimum \ 2 vs. ibi n at fixed values 
of Mi. Line colours and styles have the same meaning as in Fig. 
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Figure A3. The ratio of total i' luminosity to that contributed 
by the disk alone as a function of L g i ^ = q at <fr = 0.0 and 0.0 and 
fixed Mi . The different colours indicate the same M± values as in 
the previous two plots. 



Lgi^—Q, there is no apparent difference in best-fitting ibin 
between Mi values. 

The origin of the i' data preferring heavy NSs lies in the 
amplitude of i' light curves. In particular, since the March 08 
i'-data probes almost the full amplitude of the light curve, 
stronger constraints on input parameters are derived from 
it versus the March 10 data. Why is the i' amplitude im- 
portant? Figure [A"3l shows the ratio of total i' luminosity to 
that contributed by the disk at both <j> = 0.0 and 0.50 as 
a function of L g i ^ =0 (this is for the best fitting models at 
fixed Mi with the same values as in the previous two fig- 
ures). As L g i^—o increases, the disk contribution to the total 
light increases, but is typically less in the systems with more 
massive Mi. The evolution in disk contribution at ^ = 0.5 
is particularly strong. 

With a strongly varying disk contribution to the to- 
tal light at <f> = 0.5, the overall i' amplitude will vary 
with L g i t< j, = Q and this amplitude differs between Mi values. 
Typically, the best-fitting models at higher Mi and fixed 
Lgi^—Q have a larger i' amplitude. This is due to heavier 
Mi typically having more massive M% which are intrinsically 
brighter (this can be seen by the difference even at <f> = 0.0, 
where larger Mi systems have smaller disk contributions). 

The lower panel of figure lA"4l looks at how the i' ampli- 
tude varies with L g i ,^ = o and Mi for the best fitting models 
shown in Fig. IA3I It can be seen clearly there that larger Mi 
values produce larger i' amplitudes at fixed L g i^ =0 . The 
horizontal dashed and dotted lines show the approximate 
amplitude and lower la error bar for the March 08 i' data. 
The fact that the best-fitting models with larger Mi agree 
in i! amplitude with the data at larger L g i^ =0 values is the 
ultimate reason why the L g i^—o limits favour massive Mi. 
This is shown in the upper panel, where the minimum \ 2 for 
the three different Mi values is reproduced along with the 
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Figure A4. Lower panel: the amplitude in i' produced by the 
models in Fig. IA3I (coloured lines) versus the March 08 i' data's 
amplitude (dashed line with 1-a error bar given by dotted line). 
Upper panel: the x 2 of these models with the L g i j,_q limits cor- 
responding to a 3 and 10% uncertainty on the d value (vertical 
dashed and dotted lines, respectively). 

Lgi ^ =0 limits corresponding to the 3 and 10% d uncertain- 
ties. This also shows why if d were only known to 20%, there 
are no constraints on Mi : the best fitting low- Mi models at 
lower Lgi^ =0 are then not excluded. Thus, the lower L g i^—Q 
limit, given our current data, is critical to the constraints on 
Mi derived here. 
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